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Abstract. Let D C R'^hea, bounded domain and denote by V{D) the 
space of probability measures on D. Let 

L = • aV + 6V 

be a second order elliptic operator. Let /j, € V{D) and 5 > Q. Consider 
a Markov process X{t) in D which performs diffusion in D generated 
by the operator 5L and is stopped at the boundary, and which while 
running, jumps instantaneously, according to an exponential clock with 
spatially dependent intensity V > 0, to a new point, according to the 
distribution ^i. The Markov process is generated by the operator Ls^^^v 
defined by 

Ls^i^^v<P = SL(j) + V{ 4>dn-(j)). 

J D 

Let (p6,fi,v denote the solution to the Dirichlet problem 

L5,i^,v(t> = in D; 
4> = f on dD, 

where / is continuous. The solution has the stochastic representation 

(f>S,^,v{x) = E^f{X{TD)). 

One has that 0o.fj.v(/) = lim^^o 4>s,ti,v{x) is independent of x G D. We 
evaluate this constant in the case that n has a density in a neighborhood 
of dD. We also study the asymptotic behavior as 5 of the prin- 
cipal eigenvalue \o(5,fi,,V) for the operator Ls^fi^v, which generalizes 
previously obtained results for the case L = ^A. 



1. Introduction and Statement of Results 

Let D C he a bounded domain with C^'"-boundary {a G (0, 1]) and 
let V{D) denote the space of probabihty measures on D. Let 

(LI) L = ^V-aV + 6V 

be a second order elliptic operator. Assume that the coefficients a = {ajj}"^^ 
and b = {bi}2=i are in C^'°'{D) and that a{x) is positive definite for each 
X e D. Fix a measure n G V{D) and fix 5 > 0. Consider a Markov process 
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X{t) in D which performs diffusion in D generated by the operator 5L and is 
stopped at the boundary, and which while running, jumps instantaneously, 
according to an exponential clock with spatially dependent intensity V, to a 
new point, according to the distribution /u. That is, the probability that the 
process X{-) has not jumped by time t is given by exp(— /q^^^ V{X{s))ds), 
where td = inf{t > : X{t) D} is the first exit time from D. From its 
new position after the jump, the process repeats the above behavior indepen- 
dently of what has transpired previously. We assume that V is positive and 
in C°'{D). Denote probabilities and expectations for the process starting 
from X G D hy P^^''^ and Ex^'^ . 

Let Ls^fiy denote the operator defined by 

L5^^y4> = 5L<j) + V{l 4>dfi-(j)). 
Jd 

The operator Ls^^y generates the Markov process X{t), and consequently, 
A td)) — /o^^^ Ls^f^y4>{X{s))ds is a martingale. 
Let 4'5,fiy denote the solution to the Dirichlet problem 

Ld,^iy4> = in D; 

(1.2) 

(j) = f on dD, 

where / is continuous. It follows that (l)s,fiyiX{tATo)) is a martingale; thus 
(j){x) = Ex<t)5^^y{X{t A Td)), for all t. Letting t — )• oo gives the stochastic 
representation 

(1-3) 4>s,^.y{x) = Exf{X{TD)). 

In this paper we investigate the behavior of (ps^^y as 5 — )• 0; that is, in the 
small diffusion limit. Since L^^^ycj) = V{x){f^(l)dfi — (/'(x)), one expects 
that lims^o (j)s,fiyix) will be independent oi x £ D, and we can prove this 
trivially via the stochastic representation in ()1.3p . We wish to calculate the 
constant 



(1.4) 



4>o,fiy{f) = li'al<^)5,^iy{x), x £ D. 
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The above constant depends very strongly on the behavior of /i near the 
boundary. Here we treat the case that supp(/x) H dD ^ and that ^ has a 
density in a neighborhood of the boundary. The density may vanish on the 
boundary. Let L denote the formal adjoint of L: 

L = • aV - 6V - V • 6. 

Theorem 1. Let D C , d> 1, be a bounded domain with a C'^'" -boundary 
(a £ (0, 1]) and let fi £ V{D). Assume that V > on D. Let = {x £ 
D : dist{x,dD) < e}. 

Assume that for some e > 0, the restriction of fi to possesses a density: 
fi{dx)\£)t = fi{x)dx. Assume that for some k > 0, the following conditions 
hold. If k is even, assume that /x G C'^{D'^); if k is odd, assume that ji G 
C'^+HZ)'). Assume that 

dl^fi 



dx^ 
d^fi 



= on dD, for all < /c - 1, if k > 1; 



^ on dD, for some \(3\ = k, ifk>0. 



dxP 

Assume that V G C2'°(Z)), if k = 0,2, and that V G C^{D), if k > A is 
even; assume that V G C^~^^{D), if k is odd. 

Assume that aij,bi G C^'"(L>), if k = 0,2, and that aij,bi G C'^^^(Z)), if 
k > A is even; assume that aij,bi G C^'"(Z)), if k = 1, and that aij,bi G 
C''{D), ifk>3is odd. 

Let n denote the inward unit normal to D at dD. Let a denote Lebesgue 
measure on dD. 

If k is even, then the solution (ps.^y to (|1.2p satisfies (jl.4p with 

far, f\/ in ■ an)V^~^ L2 nda 
(1-5) ^o,.,y(/) = L-^f ^- 

If k is odd, then the solution (j)s,fj,.v to (II. 2p satisfies ()1.4p with 



Jqd fV '''^^ aV {L" ■2.' fi) ■ nda 



(1.6) 4>o,^y{f) 



JqjjV ''2^ aViL^i" h) . nda 
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In particular then, % 



if k = 0, one has 



(1.7) 



!aD ^/{n-an)^dcj 



and if k = 1, one has 



(1.8) 



</'o,M,y(/) 



IdD ^aVfi ■ nda 



IdD ^ ^aVfJ, ■ nda 



Remark 1. Note that if A: = or A; = 1, then (po,fi,vif) depends on the 
diffusion coefficient a, but not on the drift coefficient b, whereas for k > 2, 
(po^fiyif) depends on a and b. 

Remark 2. If fi has compact support, the behavior of 4'o,fj,,vif) is com- 
pletely different. In this case, (po,fi,v{f) can be studied using Wentzell- 
Freidlin action functionals. Assuming the uniqueness of the minimum of a 
certain such functional, (f>o,ny{f) will have the form /(xq) for some xq E dD. 

Define the contraction semigroup 



where Cq{D) is the space of continuous functions on D vanishing on dD. 
The infinitesimal generator of this semigroup is an extension of the oper- 
ator Ls^^y, defined on C'^{D) n {(j) ■ (p,Ls^^y(j) ^ Co{D)} with the ho- 
mogeneous Dirichlet boundary condition. The operator T^'^'^ is compact 
(see [3] where the case of constant coefficients is considered); thus, the 
resolvent operator for T^'^'^ is also compact, and consequently the spec- 
trum a{Ls^^y) of Ls^^y consists exclusively of eigenvalues. By the Krein- 
Rutman theorem, one deduces that —Ls,^y possesses a principal eigenvalue, 
Xo{d, ^,V); that is, Xo{5,n,V) is real and simple and satisfies Xo{6,fi,V) = 
inf{Re(A) : A E 0"(— L^^^^y)} [5j. It is known that A E a{—Ls^^y) if and 



r/''^'^/(x) = i5;f'^'^(/(X(t));rz5 > t), / E Co{D) 
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only if exp i-Xt) G a{Tt''^''') i2j. Thus, since HTj'^' || < 1, it fohows that 
Xo{6,fi,V) > 0. We have 

sup ||T/''^'^/|| = supPi''^'^(rz5>t); 
/eCo{5),||/||<i xGD 

thus, a standard result |6J allows us to conclude that 

lim I log sup P^'f^'^iTD >t) = -Ao(<5,/x,y). 

i-s>oo t X&D 

It is well known that this is equivalent to 

(1.9) hm ilogPi'^'^(r^ >t) = -Xo{S,fi,V), x e D. 

t— >oo t 

In the case that L = ^A, that is the case that the underlying motion is 
Brownian motion, the papers [3], [3] investigated the behavior of the prin- 
cipal eigenvalue Ao(5, /U,F) as 5 — 0. (Actually, in those papers, one finds 
the operator jLi ycj) = ^A(/) + "fV{f^ (pdfi — (j)) with 7 — )• 00.) The key 
calculations contained in Proposition [T] of this paper for the case of a gen- 
eral diffusion operator L generalize calculations in [4J for the operator ^A. 
Using the methods of [4] along with Proposition [H one obtains the following 
generalization of the results in [4]. 

Theorem 2. Let the assumptions of TheoremUl be in effect for some k > 0. 
Then the principal eigenvalue Xq{6, fi,V) of the operator —Lg^^ y behaves 
asymptotically as follows: 
i. If k is even, 

r, f- — fc+1 ~ k_ 

fe+1 , , \(ir, \/ [n ■ an)V 2 L^iida 

(1.10) hmr — Ao(<^,^,F) = ^^^^^ ^ 1^ 

a. If k is odd, 

(1.11) im<5- — Ao(5,/i,y) = ^-^ . \, ^' . 

Remark 1. Note that if /c = or /c = 1, then the leading asymptotic 
behavior of Ao(^, V) depends on the diffusion coefficient o, but not on the 
drift coefficient b, whereas for A; > 2, it depends on a and b. 
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Remark 2. We note that if ^ has compact support, then there exist con- 
stants ci,C2 > such that exp(— C25~2) < Ao((5, < exp(— ci(5~2 ), for 
small (5 > 0. This was proven in [3, 4| for the case L = ^A. The same type 
of proof works for general L. 

If /X ^ on dD, then Theorem [2] gives 



(1.12) Ao(<5, /i, V) ^ 6-2 , as 6^0. 

Theorem [2] is proven under the assumption that V > in D. This condition 
is essential. Note that if V vanishes in a sub-domain A C D, then as 
long as the process X{t) remains in A, it never jumps, and thus starting 
from a point in A, the probability that X{t) does not exit D by time t 
is greater than the probability that a (5-L-diffusion process does not exit A 
by time t; thus, in light (jl.9p and the corresponding equation for the 5L 
diffusion process, it follows that Xo{6,fj,,V) < 6Xq{L), where A^(L) > 
is the principal eigenvalue for —L in A. In particular, Xo{6,fj,,V) is on a 
smaller order than in (jl.l2p . 

Now consider the case that V > in D but V = on dD. On the one 
hand, since the process needs to not jump in order to exit D, allowing the 
jump mechanism to weaken at the boundary should help the process exit. 
Thus, if ;U ^ on dD, one might expect that \o{S, fi, V) will be on a larger 
order than 6^ . But on the other hand, if y = in D*^, then by the argument 
in the previous paragraph, Xo{6, fi, V) < A^' {L)6. Thus, we expect that if V 
vanishes identically on the boundary to high enough order, then Xq{5,^, V) 
will be on a smaller order than (5 2. 

Now let Ve = e + where V > in D and = on dD, and substitute 
Ve for V in the righthand side of (I1.12|) . If V vanishes to the first order 
on dD, then the right hand side of (jl.l2p is on the order (e2loge)~2; in 
particular, it converges to 00. This suggests that if V vanishes to first order 
on dD, then Xo{5,iJ,,V) will be on a larger order than 6^ . li V vanishes to 



8 ROSS G. PINSKY 

second order on dD, then the right hand side of (I1.12p stays bounded and 
bounded from as e — )• 0. This suggests that if V vanishes to second order, 
then Ao(5, /i, V) will be on the order 6^ , as in (jl.l2p . If V vanishes to third 
order on dD, then the right hand side of (jl.l2p goes to as e — )• 0. This 
suggests that if V vanishes to third order or higher, then Xo{6, fi, V) will be 
on a smaller order than 5 2. 

Open Question: Consider the case that fi ^ on dD, so that if V were 
strictly positive in D, then (jl.l2p would hold. Assume that V > in D and 
that V vanishes identically on dD to the order k, k > 1. At what order does 
Xo{S, /U, V) approach when 5 — )• 0? 

In section 2 we present several auxiliary results which then allow for a 
quick proof of Theorem[TJ The proof of one of the auxiliary results is deferred 
to section 3. 



2. Auxiliary Results and Proof of Theorem [T] 

In this section we present three lemmas and one proposition, from which 
the theorem will follow quickly. We begin however with a useful construction 
of the process X{-) up to its exit time from dD. On a common probability 
space with probability measure , let {y„}^Q be an independent sequence 
of diffusion processes, where each Yn{-) is a diffusion corresponding to the 
operator 6L and stopped upon reaching the boundary, and where Yq{0) = 
X ^ D, and the distribution of 1^(0) for n > 1 is /i. (We have purposely 
suppressed the dependence of on x. Note that Yn does not depend on 
X for n > 1.) Let J-'t^n = cr(l^(,s),0 < s < t) denote the filtration up 
to time t for the process Yn{-). Denote the exit time of l^(-) from D by 
TD^n- Let Jn be a stopping time for Yn{-) satisfying V^{Jn > t\^t,n) = 
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exp(— /q^^^'" V{Yrt{s))ds). Now define by induction: 
X{t) = Yo{t), for < t < Jo A TD,o] 
if Jn-i < TD,n-i, then 

n— 1 n— 1 n— 1 

X{t) = Yn{t - Jk), for -/fc < i < X] '^i' + ^D,n- 

k=0 k=0 k=0 

Let es^xi') = Px'^'^ {X {td) G •) denote the exit distribution of the process 
X{-) from D starting from X{0) = x. Let J denote the first jump time of 
X{-). Since the distribution of converges to the point mass at oo as 5 — )• 
0, we have hm^^^o ^'^(''"D.o < -^o) = 0; equivalently, limg^Q Px^'^ (td < J) = 
0. Since Px'^'^ {X{td) G -IJ < td) is independent of x E D, we conclude 
that the weak limit of e^^x (which will be shown to exist) is independent of 
X & D. Indeed, from the above considerations and the above construction 
of X{-), we have 

lime5,^.(-) = limV^{Yi{TD,i) G -^0,1 < Ji). 
(5^0 <5— >0 

(Note that the 1 appearing in four places on the righthand side above can 
be replaced by any n >2 without changing the value of the expression.) 

For each k = 0, 1,---, let eg(-) denote the probability measure on dD 
with density CqIx) given by 

eQ{x) = ( / (n • an)V~^ fida)~^ y^{n ■ an){x)V^^ {x)L^ fi{x), 
JdD 

k > 0,k even; 

el{x) = {l V~^aV{L^ n) ■ nda)~^ {x)a{x){V{L^ fi) ■ n){x), 

JdD 

k > l,k odd. 

From (jl.3p . to prove the theorem we need to prove that if /x satisfies the 
conditions of the theorem for a particular A; > 0, then es^x converges weakly 

to Cq. 

From now on we assume that fi satisfies the conditions of the theorem 
for a particular A; > 0. Fix m > 1 and let {Aj}^^^ be a partition of D 
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into m + l disjoint connected sets satisfying tlie following conditions: (i) Aj 
has a nonempty interior, for all j; (ii) Aj D dD has a nonempty interior in 
the relative topology of dD, for all j / m + 1; (iii) 65(^1 n dD) > 0; (iv) 
dist{Ajji_^.i, dD) > 0. To prove the theorem, it is enough to show that 

(2.1) lime'liAindD) = e'^(AindD). 

Let usy denote the solution to 

SLu -Vu = OmD; 

(2.2) 

u = 1 on dD. 

Let denote expectations corresponding to V^. As is well-known, usy has 
the stochastic representation 

usyix) = S\expi- / V{Yi{t)))\Yi{0) = x). 
Jo 

Lemma 1. 

(2.3) r\Yi{0) G Aj,TD,i < Ji) = / usy dfx, j = !,■■■, m + l. 
Proof. For x E D, we have 

POO 

V\td,i < Ji\Yi{0) = x)= r\Ji > t,TD,i = dt\Yi{0) = x) = 

Jo 

fOO 

s^r^iJi > t,TD,i = dt\Tt,i,Yi{o) = x) = 

/ S\eM- V{Y^{s))ds)ldt{TD,i)\Yi{0) = x) = 

Jo Jo 

£{- [ '''' V{Yi{t))dt\Yi{0) = x) = us,v{x). 
Jo 

The result follows from this. □ 

Lemma 2. There exists a constant c = c{a, b, d,V) > depending on the 
coefficients a, b of L, on V and on the dimension d such that 

distix, dD) 

us,v{x) < cexp( J ). 

c(52 



Jo 
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In particular then by LemmaUi for some c\ > 0, 
(2.4) V\Y^iO) € A^+^,To,l < Ji) < ci exp(- 

Proof. Let Vmin = ^^^xgd ^(^)- From the stochastic representation of usy, 
we have for any t > 0, 

US,v{x) < eM-tVmin) + r^TD,! < t\Yi{0) = x) < 

(2.5) 

exp(-tymin) + ^ ( max |Yi(s) - x\ > dist(3;, aZ))|yi(0) = x). 

0<s<t 

Since L has been written in divergence form, the drift of Yi is 6{b + • a). 
Let B = max^g^ \ b{x) + • a(x)| and let A = max|.„|=i max^^f){v,a{x)v). 
From [5, Theorem 2.2-ii], it fohows that 
(2.6) 

max \Yi(s) - x\ > X\Yi(0) = x) < 2dexp(- ~ ff ), for A > 6Bt. 

Letting A = dist(x,9D) and t = {(i\st{x,dD))5^'2 in the above inequahty 
and substituting the resulting estimate on the right hand side of (|2.5p . we 
obtain 

'U<5,y(a::) < exp(-(5~^dist(x, 9D)Vtnin)+2dexp(-5"^dist(x, SZ))-*'''^ ^ ^ 



2d A 

for small S > 0, from which the lemma follows. □ 

The key result for proving Theorem [1] is the following proposition, whose 
proof is postponed to section 3. 

Proposition 1. Let the assumptions of Theorem [I] be satisfied for some 
k > 0. Let j E {1, • • • , m}. If k is even, then 

k+l f 1 /" /7 7 fc+1 ~k 

lim 2 / usydfi=—^ / \/ (n ■ an)V 2 Lijjda. 
"5^0 JAj ' v2 JAjndD 

If k is odd, then 

fc+i f 1 /" fc+i ~ fc-i 

lim 5 2 / usvdfJ' = 7; / V aV[L 2 fj,) ■ ndcr. 
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Lemma 3. Let fi satisfy the assumptions in Theorem[l\for some k >0. Let 
j G {1, • • • , m} be such that e^^Aj n dD) > 0. Then 

limV^{YiiTD,i) G Aj n dD\Yi{0) G Aj,TD,i < Ji) = 1. 

Proof. Define 

U5y,j{x) = £\l9D^A,{Yi{TD))eM- [ '''' V{Yi{t)))\Yi{0) = x). 

Jo 

An argument just like that used in the proof of Lemma [T] shows that 

(2.7) V^iYiiO) G A,,TD,i < Ji,YiiTD,i) edD- Aj) = [ usyj dfi. 
An argument just hke that used in the proof of Lemma [2] shows that 

(2.8) usy^j{x) < cexp(- 



dist(x,az? - Aj), 



Let A'' > be a positive integer and let Aj' = {x ^ Aj : dist(x, dAj) < -^j. 
By Lemma [1] and (j2.7p , we have 

r^{Yi{TD,i) GdD- Aj\Yi{0) G A„TD,i < Ji) = 

(2.9) / 1 usyjdn + f 1 Usyjdfl 

j ^ 3 



Proposition [1] of course also holds with Aj replaced by Aj^ . Using the fact 

_L_ 

that usyj < usy, applying Proposition [1] with Aj and with Aj^ , and using 

(j2.8p . we obtain 

(2.10) 

limsupP^(yi(rz5,i) €dD- A,\Yi{0) G Aj,TD,i < Ji) < . 
5^-0 eQ{Ajr\dD) 

Letting A^ — )• oo completes the proof of the lemma. □ 



We can now prove (j2.ip , which will complete the proof of Theorem [TJ 
Recall that by assumption 69(^1) > 0. Thus, by Lemmas [1] and [3l it follows 
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that 
(2.11) 

^"^(^1(0) G AuYi{td) G AindD,TD,i < Ji) = (1+0(1)) / usyd^i, as 6^0. 

JAi 

By Lemmas [J and [3] and Proposition [H it follows that 
(2.12) 

V\Yi{0) G Aj,Y^{td) G AindD,TD,i < Ji) = o(5^), for j g {2,-- - ,m}. 

Using (pTT]) . (pJ2]) and Lemma [21 we have 

es{A, n 9L») = V^iYi{TD,i) G ^1 n aZ)|rz5,i < Ji) = 

V\Yi{td) £AindD,TD,i < Ji) ^ 
VS{td,i < Ji) 

(2.13) ^y^i yS^Y.iO) G .4,, yi(rg) G n aZ), Tfl,i < Ji) _ 

, as d 0. 

Ejli Xi^ ^^5,v c?/^ + o((5 2 ) 

Now from (j2.13p and Proposition [H we have 

(2.14) limesiAindD) = e'^{Air\dD). 
(5—^0 

□ 

3. Proof of Proposition [T] 

For the proof of the proposition in the case of even k, we wih need the 
following lemma. 

Lemma 4. Let n denote the unit inward normal to D at dD. One has 

lim 6^{n ■ aVusv){x) = —\/2V{x){n ■ an){x), uniformly over x G dD. 

Proof. The proof of the result in the case that L = |A was given in [3]. In 
the proof, it was shown that everything could be reduced to local consider- 
ations. In particular, it was enough to prove that the above equation holds 
pointwise under the assumption that the boundary had constant curvature. 
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We can thus make the same assumptions here, and we can also assume that 
L and V have constant coefficients. More specifically, note that L is given 
in non-divergence by 

1*^ ^2 1 da - ■ d d 

1,1 = 1 *,7 = 1 7 = 1 

(3.1) 



i,j=l j=l 



Thus for the proof we may assume that L = \ Yji,j=i dxidxj +I^»=i -^i alj' 
where a^.j and are constant. Similar to what was done in |3j, for zero 
curvature, we take Z? = {x e ii"' : < xi < 1} and consider {n ■ aVusy){0); 
for curvature i? > 0, we assume that D = Ar ^ = {x G R'^ : < |2;| < R} 
and consider {n ■ a'Vusy){x), for some x with |x| = R; and for negative 
curvature —R < 0, we assume that D = Aji^2R = {x ^ : i? < |x| < 2R} 
and consider (n • a^usy){x), for some x with |x| = R. We will consider the 
cases of zero curvature and positive curvature; the case of negative curvature 
being handled similarly to the case of positive curvature. 

We begin with the case of zero curvature. Let {cj}^^^ denote the standard 
basis vectors. Let Hi denote the hyperplane x\ = 0. The interior unit 
normal to D on dD n Hi is constant and equal to ei; that is, 

n = ei. Let y denote the projection of n-a onto Hi. Then y + {n-an)ei = 
n ■ a. Since usy = 1 on Hi, we have 
(3.2) 

(n . aVn, v-)(0) = lim ^syjtn ■ a) - usyjO) ^ usyjtn ■ a) - usyjty) 

v usy{ty + t{n- an)ei) - usy{ty) Busy 

lim = n ■ an lim — [ty + st[n ■ anjei) 

t->-o+ t t^o+ OXi 

(n • an){n ■ VuAy)(0), 



where < < t. Since usy depends only on xi and since n = ei, we 
can reduce the calculation of {nVusy){0) to a one-dimensional problem. 
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So we write usy = usy{x) with < x < 1. Now u solves the con- 
stant coefficient equation ^6ai^iu'^y + SBiu'^y — Vu = with the bound- 
ary condition n5y(0) = y(l) = 1. The quantity {n ■ Vn^ y)(0) above 
is now given by u^y(O). One can solve this explicitly and check that 
lim^^o ^"^u'^y (0) = — \J Substituting this in (I3.2p and noting that 
^1,1 = n ■ an, we obtain \\mi^^Q{n ■ aVusy){Q) = —^Jl{n ■ an)V. 

Now we turn to the case that the curvature is i? > 0. We let D = 
Ar n and consider the boundary point Rei. We need to evaluate lim5_^o("- " 

2 

aVusy){Rei). We first reduce the calculation to the calculation of the 
normal derivative, similar to (|3.2p . Note that the inward unit normal n = 
n{Rei) at Rei satisfies n = — ei. For small t > 0, let zt denote the point 
on l^l = R which is closest to Rei + tn ■ a. Define the vector wt by zt + 
wt = Rei + tn ■ a. (Note that Rei, zt and wt take on the roles played 
by 0, ty and t{n ■ an)ei respectively in the case of zero curvature.) Of 
course limj_^o+ — R^i\ = 0. Since the curvature is positive, we have 
\wt\ < (n • an)t; however limj_i>o+ = n ■ an. Note also that the direction 
1^ of Wt approaches the direction of n as t — )• O'^ . Thus, since ugy = 1 on 
|x| = R, we have 

f V7 \^r> \ V usy{Rei + tn ■ a) - usy{Rei) 

[n ■ a\/usv)\R^i) = hm = 

t^o+ t 

(3.3) ^.^^ Usyjzt + Wt) - Usy(zt) _ \wt\ Usyjzt + Wt) - Usyjzt) _ 

t^0+ t t-s>0+ t \wt\ 

(n • an){n ■ Vu^ y)[Rei). 

We now consider (n • Vusy){Rei). Let {r,0) with 9 G S'^~^ denote po- 
lar coordinates. We rewrite the constant coefficient operator L in polar 
form. Of course now the operator will no longer have constant coeffi- 
cients; however by the localization mentioned above, we may consider in- 
stead the constant coefficient operator obtained by evaluating the coefficients 
at Rei. Call the resulting operator L. We have L = |(n • an)-^ + B-^ + 
terms involving differentiation with respect to 6 and maybe also r, where 
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i3 is a certain constant whose form is irrelevant for our purposes. Now 
usy solves 6Lusy — Vugy = for ^ < r < R, and ugy = 1 at r = 
and r = R. It follows that usy is a function of r alone. Thus usy satis- 
fies the one-dimensional equation ^5{n ■ an)u'^y + SBu'^y — Vu^y = for 
^ < r < R and u{^) = u{R) = 1, and (n • 'Vusy){Rei) becomes —u'^yiR). 
We have thus reduced the problem to the previous case of zero curvature, 
and conclude that lim5_5.o("' ■ a''^usy){Rei) = —y^2{n ■ an)V. 

□ 

Proof of Proposition 1. Let /io(") be an arbitrary probability measure on D 
which has a density ^q{x) which satisfies the same smoothness assumptions 
in D that the density /i satisfies in Z)'^, and which satisfies the same vanishing 
conditions on dD that the density /x satisfies there. An easy argument then 
shows that to prove the proposition, it suffices to prove it with Aj replaced 
by l), d/x replaced by d/xo and /x(x) replaced by ^io{x). We will first prove the 
proposition for the case fc = 1, which is easier than the case /c = 0. We then 
show how to go from the case A; = 1 to the case A; = 3, from which it will be 
clear how to proceed for odd k. After that we will prove the proposition for 
A; = and then we show how to go from the case /c = to the case k = 2, 
from which it will be clear how to proceed for even k. 

In light of the above paragraph, we consider f^usyfj-odx. Since k = 1, 
/io vanishes on dD, but V/xq does not vanish identically on dD. Using (|2.2p 
and the fact that /xq vanishes on dD, and recalling that n denotes the inward 
unit normal, integration by parts gives 

/ usvfJ'odx= / Lusv^dx = 
(3.4) ^Jd V 

/ usyL^dx + - a\7{^)-nda. 
J D V I Jqq V 

(Note that by assumption, /io and V are C^-functions so there is no problem 
with the integration by parts.) By Lemma [21 usy converges to boundedly 
pointwise on D. Also, since /xq vanishes on dD, we have V(^)-n = yV/XQ-n 
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on dD. Thus, letting (5 — )• in (I3.4p . we obtain 

We now turn to the case k = 3. In the case A; = 3, and all its derivatives 
up to order 2 vanish on dD; in particular, the last term on the right hand 
side of p.4p is 0. Thus, using (|2.2p again, integrating by parts and using 
the fact that the second order derivatives of fiQ vanish on dD, we have from 

Jd J d V J £, V V 

I usyL^L^dx + / aV{^L^) ■ nda. 
JD V V Jqq V V 



(3.5) 



(Note that by assumption, hq and V are C^-functions and a,j- and 5j are 
C^-functions, so there is no problem with the integration by parts.) Using 
Lemma [2] again and the fact that and all its derivatives up to order 2 
vanish on dD, we obtain 

limS^^ / usvfJ-odx= / V'^aV {L^q) ■ nda. 

The same technique is used repeatedly to handle larger values of odd k, the 
smoothness requirements in the statement of Theorem [1] being the smooth- 
ness required to implement the integration by parts. 

Now we turn to the case k = 0. Let w solve the equation 

~ w 

L- = in D; 

(3.6) ^ 

w = fj-o on dD. 

Note that by the smoothness assumptions on a, b, V, fio, it follows that w is 
the solution to an elliptic equation with C"^(Z))-coefficients and continuous 
boundary data. Thus, w G C^'°'{D) D C{D). 
We will show below that 

(3.7) lim(^~2 / Us vifJ-o ~ '^)dx = 0. 
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Thus, it is enough to show that 

(3.8) hm(5^2 / usywdx = —= f ^ in ■ an)-^^da. 



JM)_ 



Using (j2.2p and (j3.6p . and integrating by parts, we have 

(3.9) §-2 f usvwdx = 62 [ {Lus v)j7dx = f f^aVusv ■ nda, 
Jd J d ^ JdD ^ 



where we have used the fact that 

If w /"^, f ~ w 

- 1 aV( — ) • nd(J — I —b ■ nda = / L—dx = 

2 JdD y JdD y Jd V 

by (|3.6p . (Actuahy, since w is not necessarily up to the boundary, in 
the above integrals one should replace D hy D — D'' and dD by d{D — D'') 
and then let e — )• 0.) Letting 5 — )• in (|3.9p . and using Lemma HI we obtain 

(IHSl). 

It remains to prove (|3.7p . By Lemma [21 we have 
(3.10) 

|(5~2 / usy{lJLQ — 'w)dx\ < snY>{lJLQ{x) + w{x))\D\5~'iceyiY>{ -). 

Jd-d^ ' x&D c52 

We also have 

(3.11) |(5^2 / usy{fio — w)dx\ < sup \fio{x) — w{x)\{S^^ / usydx). 
Jd' ' xeD^ Jd 

Now (13. 8p holds for every fiQ in a wide class; in particular, it holds for fiQ 
which are uniformly positive on D. In such a case, it follows by the maximum 
principal that w is uniformly positive on D. (The principal eigenvalue for L 
coincides with that of L, and is consequently negative. Thus the generalized 
maximum principal holds: Lv = in D and w > on dD guarantees that 
V > on Z). Apply this with v = y.) By considering (]3.8p with such a 
uniformly positive w, it follows that 5~2 usydx is bounded as (5 — 0. 
Using this, the proof of (13. 7[) now follows from (13.10p . (I3.11j) and the fact 
that lim^^o sup3^.g£)e |^o(^) ~ w{x)\ = 0. 
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We now turn to the case k = 2. Since and all its derivatives up to 
order one vanish on dD, we can write (13.41) as 



(3.14) lim I usy{L^ - w)nodx = 0. 



(3.12) (5 2 / usyfiodx = 52/ usyL^dx. 

Jd ' Jd ' ^ 

As with the case A; = 0, we define an auxiliary function w. Let w solve the 
equation 

L— = in D; 

(3.13) ^ 

w = L— on oU. 

(By assumption, ^0 and its first order partial derivatives vanish on 9D, but 
not all of its second order partial derivatives vanish on dD. It then follows 
from the maximum principal that > on dD.) The same argument 
used to show (|3.7p shows that 

2 „ _ 

In light of (j3.12p and ()3.14p . it is enough to prove that 

(3.15) lim(5^2 / usvwdx = —= i \/ (n ■ an)V~^LuQda. 

<5^o Jo ' V2JdD 

Using (j2.2p . integrating by parts and using (j3.13p . we have 

(3.16) 

(5~2 f usywdx = 52 f {Lusy)^dx = f ^{L^)aVusy ■ nda, 
Jd J d V z Jqj^ V V 

where we have used the fact that 

If f w f 

- I a\/[—)-nda— / —b-nda= I L—dx = 

2 JdD V Jqd V Jd V 

by ()3.13p . Since /xq and all its first order partial derivatives vanish on dD, 

we have = yLfiQ on dD. Using this and LemmaUl and letting 5 — )• in 

(j3.16p . we obtain (j3.15p . The same technique is used repeatedly to handle 

larger values of even k, the smoothness requirements in the statement of 

Theorem [1] being the smoothness required to implement the integration by 

parts. □ 
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